High fidelity simulations of ion trajectories in miniature ion traps 
using the boundary-element method 
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In this paper we present numerical modeling results for endcap and linear ion traps, used for 
experiments at the National Physical Laboratory in the UK and Innsbruck University respectively. 
The secular frequencies for S8 Sr + and 40 Ca + ions were calculated from ion trajectories, simulated 
using boundary-element and finite-difference numerical methods. The results were compared against 
experimental measurements. Both numerical methods showed high accuracy with boundary-element 
method being more accurate. Such simulations can be useful tools for designing new traps and trap 
arrays. They can also be used for obtaining precise trapping parameters for desired ion control when 
no analytical approach is possible as well as for investigating the ion heating rates due to thermal 
electronic noise. 

PACS numbers: 03.67.Lx, 32.80.Pj 



I. INTRODUCTION 

The manipulation of laser-cooled, trapped ions within 
radio-frequency (RF) traps has been widely studied for 
atomic optical frequency standards 0, 0- In the last 
decade, ion traps have also been used for quantum in- 
formation processing. This began with the first general 
method for the implementation of quantum-logic gates 
with ion traps and hardware realization of the C-NOT 
gate Q . Then the multiplexing trap scheme [f| was pro- 
posed to improve the scalability, which is an important 
issue for quantum computing. This was followed by an 
architecture for a large-scale ion-trap quantum computer 
6] and a dual linear trap array for transferring ions from 
one trap to another to allow sequential quantum-logic 
operations to be performed \J\. 

More recently, research has turned towards the minia- 
turization of ion traps and investigation of different elec- 
trode geometries for microfabrication of trap arrays 0,0. 
Typically, the trap electrodes are suggested to be planar 
since they are much easier to fabricate than circular and 
hyperbolic geometries, especially at such a small sizes. 
These micro-trap arrays would allow the construction of 
more complicated devices for ion trapping at microscopic 
level, which should satisfy some of the main requirements 
for scalable quantum computing. Such miniaturization 
techniques can also be used for the size reduction of 
miniature quadrupole mass filters 0] ■ Quadrupole mass 
filters are mass analyzers like ion traps and they can also 
be constructed to form arrays. 

Analytical and numerical potential modeling is crucial 
for ion trap design since static and ponderomotive poten- 
tials are responsible for ion oscillations within the trap. 
Numerical modeling of ion trap electrostatics has already 
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been done by using different approaches, which include 
finite-difference method (FDM) jll| and finite-element 
method (FEM) @. Both FDM and FEM can produce in- 
accurate potentials and fields, particularly near the edges 
of electrodes where the results can be highly inaccurate. 
Such inaccuracy near electrode edges could have an effect 
on mass spectra in ion trap mass spectrometers where 
ions often approach electrode edges. High accuracy can 
also be an important factor for performance predictions 
in complicated electrode structures like trap arrays. Even 
if small inaccuracies are made in calculating suitable elec- 
trode dimensions and distances, the errors can accumu- 
late for multi-zone trap arrays [lj • This could later have 
an effect on ion transfer between different regions and 
ion motional heating. For these reasons, the boundary- 
element method (BEM) has also been considered for nu- 
merical modeling of quantum computing ion traps |13| . 
Unlike FEM and FDM, which use the whole electrode 
volumes to define grid points, BEM uses only the surface 
of electrode volumes since there is no need to consider 
the effects of charges below the surface. This enables 
faster computation and higher accuracy even with small 
number of electrode segments. The direct comparison of 
BEM, FDM and FEM has been demonstrated in Cubric 
et al. 0| for an ideal spherical analyzer and double cylin- 
der lens using different benchmark tests. These tests in- 
cluded simulations with commercially available programs 
CPO (BEM) [HI and S1M10N (FDM) 0. All the re- 
sults showed that BEM had much smaller error levels 
than FDM and FEM for modeling potentials, fields and 
particle trajectories. 

In Sec. II and Sec. Ill, the ion secular frequencies were 
calculated from numerically obtained ion trajectories for 
both an endcap and circular linear trap. Endcap traps 
are used for atomic clocks and optical frequency stan- 
dards at NPL [17J. Circular linear traps were used for 
quantum computing experiments at Innsbruck 18]. Sim- 
ulation results were produced from CPO and SIMION 
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FIG. 1: (Color online) Schematic diagram of NPL endcap 
trap. 



and compared with experimentally measured secular fre- 
quencies 

[3 S3, E3- CP0 produced closer results to 
measurements than SIMION for both traps. The heat- 
ing rate due to Johnson noise was also calculated for the 
Innsbruck linear trap. 




FIG. 2: (Color online) Equipotential contours in zx(zy) plane 
for the NPL endcap trap. 

and y with r (= \J x 1 + y 2 ). Thus, the trap potential at 
any point in a conventional quadrupole ion trap is given 
by: 

Hr, z) = ^ ~ 2 f 9 + 2 2Z °V + Vcos(fit)), (1) 

where U and V are the DC voltage and zero-to-peak AC 
amplitude applied to the ring electrode, f2 is the angular 
frequency equal to 2nf, where / is the frequency of the 
RF field, ro is the smallest distance from the trap center 
to the ring electrode and zq is the smallest distance from 
the trap center to the endcap electrode. For an ideal 
quadrupole field the ratio between ro and zq is given by: 



II. ENDCAP TRAP 

The endcap trap, proposed by Schrama et al. |23, is a 
different geometric variant of a conventional quadrupole 
Paul trap and is a non-linear trap. Fig. 1 shows a 
model of the endcap trap used at NPL for atomic clocks 
and frequency standards experiments with 88 Sr + and 
171 Yb + [23. A conventional quadrupole ion trap consists 
of one ring electrode (the hyperbolic cylinder) to trap 
ions in radial direction and two endcap electrodes (hy- 
perbolic plates) for capturing ions in the axial direction. 
To produce an ideal quadrupole field, an RF potential is 
applied to the ring electrode, while the endcap electrodes 
are grounded. As it can be seen from Fig. 1, the NPL 
endcap trap has two inner endcap electrodes and two 
outer endcap electrodes concentric with inner ones. The 
outer endcaps are moved back to allow sufficient space for 
laser access. Equal RF voltages are applied to the inner 
endcaps, while small DC voltages can be applied to the 
outer endcaps, which are normally grounded. In effect, 
the inner endcaps confine ions in the axial (z) direction, 
while the outer endcaps prevent ions from escaping in the 
radial (x,y) direction. The equipotential lines in zx/zy 
planes are shown in Fig. 2. 

The mathematical theory of an endcap trap can be de- 
scribed through the analysis of a conventional quadrupole 
ion trap. Since the electric field in quadrupole ion traps 
is rotationally symmetric, we can represent radial axes x 



Tr> — 2z n . 



(2) 



The equations of motion for an ion at mass m and 
charge e are given by: 



+ 7^s( u + Vcos(Qt))r = 0, 
z 



- z 



dt 2 



■(U + Vcos(nt))z = 0. 



(3) 



Ion stability parameters a u and q u are obtained by solv- 
ing the Mathieu equation: 



— + (a u - 2q u cos(2£))u = 0, 



(4) 



where u can be either x, y or z and £ = (cot)/2. The 
resulting expressions for a u and q u are: 
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(5) 
(6) 



where e is the 'efficiency' of the trap. A conventional 
quadrupole ion trap produces an ideal quadrupole field 
and it has s = 1. It was shown experimentally that the 
NPL endcap trap has e = 0.63 [2!^. 
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TABLE I: Experimental and numerical results for secular frequencies of a 88 Sr+ ion, trapped within the NPL endcap trap. 



Secular frequencies 


Description 


1st set 


2nd set 


3rd set 


4th set 


5th set 






199 V RMS 


221 V RMS 


245 V RMS 


274 V RMS 


304 V RMS 






15.955 MHz 


15.948 MHz 


15.936 MHz 


15.925 MHz 


15.91 MHz 






Z.lZ V L)Vj 


Z.00 V D\y 


o.ol V LJb 


o qo \r tic 


Z.oo V 1JO 


LU Xi y/2lT 


Experiment 


1.395 MHz 


1.590 MHz 


1.800 MHz 


1.980 MHz 


2.230 MHz 


uj z /2tx 


Experiment 


2.985 MHz 


3.360 MHz 


3.795 MHz 


4.340 MHz 


5.070 MHz 


co Xl y/2ir 


BEM 


1.403 MHz 


1.596 MHz 


1.789 MHz 


1.980 MHz 


2.227 MHz 


uj z /2n 


BEM 


2.939 MHz 


3.265 MHz 


3.767 MHz 


4.281 MHz 


4.960 MHz 


u x ,y /2n 


FDM 


1.441 MHz 


1.606 MHz 


1.791 MHz 


1.988 MHz 


2.213 MHz 


u! z /2n 


FDM 


2.879 MHz 


3.247 MHz 


3.668 MHz 


4.261 MHz 


4.946 MHz 



An ion is stable within the ion trap if it has a stable 
trajectory in both radial and axial directions. Another 
important trapping parameter is f3 u , which depends on 
a u and q u and < /3 U < 1 must hold. In order to obtain 
the exact value of (3 U , a continued fraction in terms of a u 
and q u must be used. A simpler expression for (3 U is the 
Dehmelt approximation given by |24| : 



Pu = [a u + (q 2 u /2)] 



1/2 



(7) 



which is only valid for q XtV < 0.2 and q z < 0.4. Because of 
its trap parameters (e.g. RF voltages), the NPL endcap 
has larger values of q u and the Dehmelt approximation 
cannot be used. Therefore, the fourth order approxima- 
tion for j3 u should be used |25j : 



(5a u + 7)q* 



2(a„ - l) 2 



32(a u -l) 3 K-4) 



1/2 



Ion motion in the RF field consists of secular mo- 
tion (slow oscillations) and micromotion (fast oscilla- 
tions). In an endcap trap, the ion experiences micro- 
motion in all three directions. Micromotion can cause 
adverse effects for laser-cooled ions, such as significant 
second-order Doppler shifts when high-accuracy is inves- 
tigated and limited confinement time in the absence of 
cooling. This is due to the increase of ion motional heat- 
ing. In experiments, micromotion is normally minimized 
with the fluorescence modulation technique [26j and com- 
pensation electrodes, which are used to move the ion to- 
wards the trap center where the energy of micromotion 
is the lowest. The NPL endcap trap has two compensa- 
tion electrodes orthogonal to each other and to the trap 
electrodes, which reduce micromotion in radial direction, 
while small DC voltages can be applied to outer endcap 
electrodes for the reduction of micromotion in the ax- 
ial direction. The ion motional frequencies are normally 
called secular frequencies, since micromotion is very small 
compared to secular motion and its influence can be ne- 
glected for high frequencies. The expression for angular 
secular frequencies is an algebraic progression and it is 
given by |24j : 



< n < oo, 



(9) 



where n represents the frequency order. The angular 
secular frequency at n = is: 

u u = —, (10) 

and it is called the fundamental frequency, having the 
lowest value of all orders and the highest power spectrum. 
The expressions for ui x ,y and ui z can be obtained by using 
the approximation for f3 u from Eq. (8) and placing it into 
Eq. (10). The inner endcap DC voltage U is always set 
to zero, so that a u = for the NPL endcap trap and only 
q u is necessary for calculations of (5 U . 

To model ion traps accurately, all the trap parameters 
need to be explicitly specified. This includes both trap 
dimensions and driving voltages. The NPL endcap trap 
consists of two inner and two outer endcap electrodes 
made from Tantalum, which have an Alumina insulation 
spacer between them. The inner encaps have 0.5 mm 
diameter and a length much larger than their diameter 
(approx. 16 mm). The outer endcaps have 1 mm inner 
diameter and 2 mm outer diameter. The inner endcaps 
are separated from each other by 0.56 mm, which is equal 
to 2zo- The outer endcaps are separated by 1 mm and 
angled at 45 degrees with respect to the z-axis, which is 
also called the trap axis. In the simulations, the trap was 
driven with different voltages and at different frequencies 
corresponding to a few sets of experiments performed at 
NPL. A 88 Sr + ion was injected at the trap center at 0.05 
eV kinetic energy with equal initial velocities in all direc- 
tions and allowed to oscillate for 1 ms. 

Table I shows the numerical and experimental values 
for secular frequencies of a 88 Sr + ion for five different sets 
of experiments at NPL. In all the experiments a single 
laser-cooled ion was confined near the trap center. Each 
set contains different trapping parameters. These include 
the RF voltages and driving frequencies applied to the in- 
ner endcaps and small DC voltages applied to one outer 
endcap, while the other one remains grounded. The in- 
ner endcap DC voltages are equal to zero. The secular 
frequencies for an endcap trap cannot be obtained ana- 
lytically like for a conventional quadrupole ion trap using 
Eq. (9). This is because an endcap trap does not have 
an ideal quadrupole field and e = 1. Thus, the value of 
its efficiency can only be estimated experimentally or nu- 
merically. The secular frequencies were calculated from 
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FIG. 3: (Color online) Schematic diagram of Innsbruck linear 
trap with circular electrodes. 

ion trajectories produced by CPO (BEM) and SIMION 
(FDM), and compared with exact experimental measure- 
ments taken from |l9j| . The measurement technique used 
at NPL is described in [2]} , where radial and axial secular 
frequencies can be seen from the experimental sideband 
spectrum. 

Before modeling the NPL endcap trap, BEM and FDM 
simulations were performed on a conventional quadrupole 
ion trap to compare their results against theory (see Eq. 
(1)). These simulations showed that BEM was in av- 
erage 1% more accurate than FDM for basic potential 
and field calculations near the trap center. Table I also 
shows that on average BEM produced closer secular fre- 
quencies to the measurements for the NPL endcap trap 
than the ones FDM generated. Obviously this is a basic 
comparison and it would require more experimental sets 
and simulations to provide a definitive picture. It can be 
seen from the table that the radial frequency is smaller 
than the axial frequency, which is always the case for an 
endcap trap since the RF field is stronger in the axial 
direction due to the RF driven inner endcaps. 

From the secular frequencies, the DC component of 
the electric quadrupole field gradient can be calculated, 
which can be used to determine the quadrupole moment 
of an ion ^t|- This can be useful when designing traps 
for atomic clocks and frequency standards with different 
ions. 



III. LINEAR TRAP 

A conventional linear Paul trap consists of a 
hyperbolic-rod quadrupole mass filter with prefilter and 
postfilter of the same shapes, which represent endcap 
electrodes. The main filter traps ions radially with RF 
voltages, while the prefilter and postfilter confine them 
axially with DC voltages. Because of the difficulty in 
the manufacture of hyperbolic electrodes, other geomet- 
ric shapes are used for linear traps, such as circular-rod 
[ls|| . blade-shaped |2^|, rectilinear [2^ and planar 0, 




FIG. 4: (Color online) Equipotential contours in zx/zy plane 
for the Innsbruck circular-rod linear trap. 



which have all found applications in different areas. 

This section will investigate only linear traps that 
are used for quantum computation. Fig. 3 shows the 
circular-rod linear trap, which represents the 'old' trap 
model used for quantum computing with 40 Ca + ions at 
Innsbruck University. It has four circular rods with two 
ring-shaped electrodes placed around the ends of the 
quadrupole rods. Ions are radially trapped with RF volt- 
ages applied to the quadrupole rods, while axially they 
are confined by DC voltages applied to the endcap rings. 
In linear traps used for quantum computing one pair of 
the two diagonally opposing RF electrodes is grounded, 
while in linear trap mass spectrometers, the two pairs of 
diagonal rods have equal voltage values, but opposite in 
sign. The ratio between the rod thickness and diagonal 
rod distance is lower for quantum computing traps than 
for mass spectrometry traps. This is because quantum 
computing traps do not require ion filtering and they 
need larger space between the electrodes to allow access 
for laser beams. Also there is no DC voltage applied to 
the RF electrodes, since again ions do not need to be 
filtered, but only to remain trapped. 

Linear trap theory is based on quadrupole mass filter 
theory. The electric fields inside a linear trap have com- 
ponents along the trap axis (z direction) and radial axes 
(x and y directions). Fig. 4 shows the equipotential lines 
in zx/zy planes for the Innsbruck linear trap. The trap 
potential at any point in the radial direction is given by: 



2 2 

<f>( X ,y) = ±^y-(U+VcOs(nt)), (11) 
r 

where U and V are the DC voltage and zero-to-peak AC 
amplitude applied to the quadrupole rods, O is the angu- 
lar frequency equal to 27r/, where / is the frequency of 
the applied RF field, and ro is the smallest distance from 
the quadrupole center to the surface of the electrodes. 
Therefore, the equations of motion for an ion at mass m 
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and charge e can be derived: 



fx 
dt 2 

fy 
dt 2 



r-^-^(U + Vco S (Qt))y = 0, 



d 2 ^ 
dt 2 



0. 



(12) 



Ion stability parameters a u and q u are obtained by solv- 
ing the Mathieu equation given in Eq. (4). The resulting 
expressions for a u and q u are: 



1 AenUo 

2 z mVl 2 ZQ 2 ' 



1x 



-Qv 



mfl 2 ro 2 



fe = 0, 



(13) 
(14) 



where k is the 'geometric factor' |2y| which must be es- 
timated, Uo is the DC voltage applied to both endcap 
rings and z is the smallest distance from the trap center 
to the endcap electrode. In practice the geometric factor 
describes the intensity of the static DC field in the axial 
direction. Its value depends on the change of trap dimen- 
sions and the static DC voltages applied to the endcap 
electrodes. 

In theory, an RF field in linear traps has no compo- 
nents in the axial direction. However, when a whole 3D 
structure is modeled numerically, a small presence of RF 
field can be seen in the axial direction. This means that 
micromotion also exists in the axial direction, but it is 
much smaller than in radial. The expressions for radial 
and secular frequencies can be obtained by using the ap- 
proximation for f3 u from Eq. (7) and substituting it into 
Eq. (10). Since \a x ,y\ « \lx,y\ for the Innsbruck linear 
trap, the formula for the radial angular frequency can be 
approximated by: 



2eV 



(15) 



Since q z — 0, the expression for the axial angular fre- 
quency is given by: 



mzo 2 



(16) 



The Innsbruck linear trap consists of four stainless steel 
cylindrical electrodes, forming a quadrupole, and two 
ring-shaped electrodes at the ends. The electrodes are 
isolated by MACOR spacers. The quadrupole rods have 
0.6 mm diameters with their diagonal separation (2tq) 
equal to 2.4 mm. The endcap rings are 6 mm in diame- 
ter and the distance between them (2zq) is 10 mm. In the 
simulation, the trap was driven with 1000 V peak at 18 
MHz applied to one pair of RF electrodes with other pair 
grounded and 2000 V applied to the endcap DC rings. A 
40 Ca + ion was injected at the trap center at 1 eV kinetic 
energy with equal initial velocities in all directions and 
allowed to oscillate for 1 ms. 




a) 
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FIG. 5: BEM generated power spectra for secular frequencies 
of a 40 Ca+ ion, trapped within the Innsbruck linear trap. The 
trap driven with 1000 V AC peak at 18 MHz. 



As for an endcap trap, the secular frequencies for a 
linear trap cannot be obtained analytically because of 
the geometric factor involved. Fig. 5 shows the plots of 
40 Ca + ion motional frequencies generated by BEM for 
given Innsbruck linear trap parameters. The power spec- 
trum plot for axial frequencies has much less noise than 
the radial plot. This is because of very small presence of 
an RF field in axial direction. The experimental measure- 
ments gave uj x . y /2ir = 1.400 MHz and u> z /2ir = 700 kHz, 
which can be clearly seen from the full sideband spectrum 
given in [H,|2l|]. BEM produced lj X}V /2tt = 1.396 MHz 
and Ld z /2-K — 702 kHz, while FDM produced ui Xty /2ir = 
1.507 MHz and lu z /2it = 696 kHz. The BEM results are 
again closer to the measurements than FDM results, es- 
pecially for the radial frequency where the BEM showed 
significantly higher accuracy. It can be seen that the ra- 
dial frequency is larger than the axial, which is normally 
the case in linear traps. 

The numerical prediction of the secular frequencies can 
be used for estimating the heating rates when designing 
new traps. Ion motional heating is an important issue 
for quantum computation. It can lead to decoherence 
of superposition states and increase the ion separation 
time in multiplexed traps, which would limit the speed 
of quantum- logic operations 0, Q . The main causes of 
heating of ion motional modes include Johnson noise that 
can happen due to the resistance of the trap electrodes or 
external circuits, and fluctuating patch potentials, which 
are influenced by the noise from microscopic electrode 
regions. The heating rate due to resistance of trap elec- 
trodes is given by |30| : 



hw u n n 



e 2 kTR 



(17) 



6 



where n u is the average vibrational quantum number of 
an ion for a given direction, T is the trap operation tem- 
perature (usually T = 300 K), z is the distance of the 
ion from the conductive electrode and R is the electrode 
resistance, which can be calculated using resistivity of 
electrode material (p = 7.5 x 10~ 6 ficm for stainless steel). 
Thus, the Innsbruck trap heating rate for radial modes 
due to Johnson noise is approximately 1 quantum/670 
ms (using the BEM generated secular frequency). This is 
obviously an underestimate for the heating rate because 
of the other influential factors previously mentioned. Ion 
motional heating can especially be increased due to fluc- 
tuating potentials on electrode surfaces. However, it has 
been demonstrated that smooth and pure electrode sur- 
faces can significantly minimize the heating rate Q. If 
electrodes were made smooth enough, then the calcula- 
tion of Johnson noise heating (inevitably present) will 
give a better estimate for the actual ion heating. 

IV. CONCLUSIONS 

This paper has shown simulation results of ion sec- 
ular frequencies for two different ion traps using the 
boundary-element method and finite-difference method. 
The results were compared with experimental measure- 



ments in each case and the boundary-element method 
proved to be more accurate than the finite-difference. 
We suggest that the boundary-element method should 
be used for accurate modeling of ion-trap electrostatics 
and particle trajectories. This method should especially 
be useful when used in the design of miniature trap ar- 
rays that could potentially be used for scalable quantum 
computers. Such high fidelity simulations would help to 
prevent design errors, which could accumulate for com- 
plicated multi-zone linear traps. Ion motional heating 
is also of interest and the heating rate due to Johnson 
noise can be calculated from numerically obtained sec- 
ular frequencies. The relative permittivity of the insu- 
lating materials can be included in future simulations to 
check whether it has effects on ion secular frequencies 
and motional heating. 

For the purpose of quantum computation, a logical 
next step could be the numerical (using space-charge sup- 
port) and analytical modeling of two or more ions oscil- 
lating simultaneously in a trap. Such modeling would 
include all the mutual interactions and quantum effects 
that ions experience during their oscillations enabling the 
prediction of superposition of motional states. This could 
lead to an ion-trap C-NOT gate simulator, which might 
be further expanded to simulate quantum algorithms. 
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